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PREFACE 


A method Is presented through which estimates may be made of the 
effect of an ebaerbing, emitting » end acetterlng etmoaphere upon remote 
sensing of surface areas having non-unlfottH' Intensity. Results, presented 
In terms of contrasts and contrast transmittances, may be employed to 
predict or correct remotely sensed visible or Infrared Images. 

This work was accomplished under the technical monltorship of 
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NCmiCLAlURE 


A - Area 

C(0) - Ground level contrast 
C(h) - Contrast at altitude h 
CT - Contrast Transmittance . 
e - Hemispherical emissive power 
% - Black-body hemispherical emissive power 
h - Sensor altitude 

1 - Intensity 

Igff - Effective point-source intensity 
1 q - Ground level intensity 

R - IlLstance fron sensor line-of-sight ground intercept to effective 
point source 

2 - Altitude 

- Scattering and absorption coefficient, respectively 
TgjTg “ Optical depth, based on scattering and absorption, respectively 
t - Inclination of sensor line-of-sight with vertical 
♦ - Rotational angle locating effective point source (Figure 1) 


NCyCMCLATURE FOR APPENDIX I 


^1* ^3 “ limits on integral 

bj^, b 2 , b 3 - Upper limits on integral 

^1* ^2* ^3 “ Constats :n matrix fomulaticai (Equation [5]) 

I - Function of variables x, y, z 
I - Value of integral 

Ini “ Trapezoidal approximation of I , using ni inter'/als 
n - Number of intervals in approximate integration 

N - n + 1 

V - Volume of integration 

Wi - ith weight function in approximate integration 
x,y ,2 - Independent variables of integral 

- Values of x. v, z, respectively, at which derivatives are evaluated 
(Equations []j, [4]) 
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ATMOSPHERIC EFFECTS ON PEMOTE SENSING OF 
NON-UNIFORM TEMPEPATURE SOURCES 
FINAL REPORT 

iwrRODucrioN 

This investigation considers the effects of atmospheric participation 
ipon the remote sensing of the earth’s surface when that surface eiduhits 
a ground- level intensity vhich varies frcm point to point. This ground 
level intensity may be produced by emission at the surface, reflection, 
or both. Since the emitted energy is dependent upon the surface tenper- 
ature, it is possible to infer the surface tenperature by means of remote 
measurements, when an estijnate is made of atmospheric participation and 
other effects, to be discussed. Of particular interest in this . investigation 
is the effect of "disturbing” emitting areas, not in the sensor line-of- 
si^t, but whose energy is scattered into the sensor line-of-sight, 
pncoducing an . interfering signal at the sensor. It is the purpose of this 
investigation to provide methods of predicting remotely sensed contrast 
between two surface areas having unequal intensities, and compare that 
oontrast to the contrast at the ground between the two areas, giving a 
’'contrast transmittance". -Knowing the contrast transmittance, it is possible 
to interpret remotely sensed visible or infrared images with respect to 
resolution, loss of detail^ and other important aspects of pattern recognition. 
If prt.limlnary atmospheric measurements ir;dicate, for example, that an 
unacceptable loss of contrast would occur, a sensing run might be delayed 
pending more favorable conditions. 

CONTRIgjriNG FACTORS IN REMOTE PATTERN SENSING 

Several factors canplicate remote pattern recognition. If it is 
assumed that the upwelling radiant cner^ is due solely to tlia surface 



temperature, a larger change in energy received by the sensor will occur 
per degree of change in temperature for a warm surface than for a cool 
surface. Bastuscheck (Reference 1) has discussed this effect. The 
radiant energy received from a highly reflecting surface area, however, is 
not indicative of its tenperature. For this reason., it is preferable 
to think in terms of upwelling intensity or radiancy, with the under- 
standing that this energy my result from emission, reflection, or both. 

For the remote sensing of surface temperatures, ground- level correlations 
must be made between surface tenperat <re and total upwelling energy, taking 
into account specific emissive and reflective properties of the surface. 

Ihus the significant surface characteristic for the purposes of this 
investigation is taken to be upwelling intensity. 

Of great importance is the wavelength or frequency interval in which 
the sensing takes place. Not only does the rate of vpwelling ener^ from 
a given surface area vary greatly with wavelength, but scattering, emitting 
and absorbing effects of the atmosphere, which attenuate the signal from 
the surface, vary dramatically with wavelength. In addition, the attenuating 
effects of the atmosphere vary with location and time of year. Anding 
(Reference 2) and Chang (Reference 3) have performed detailed studies of 
atmospheric attenuation, and the 4-D atmospheric model of Chang is 
incorporated into the final resu].ts of this investigation. 

Also of great importance in this investigation is the geometry of the 
sensor line-of-sight with respect to the vertical, and with respect to 
"disturbing" emitters. 
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SCATTERING ATMOSPHERE . RESULTS FOR A POINT SOURCE 

This report is an extension of the investigation of Reference ( 4 ) . 

For continuity some of the results of that report will be included here. • 

In the previous report the purely scattering behavior of the atmosphere 

was studied wi1i\ respect to a single "disturbing” point source of kncwn 

intensity erg/sec-m^-steradian-m> not lying in the sensor line-of-sight. 

The sensor is inclined at an angle 9 frcm the vertical and an angle 

with respect to a line from the line-of-sight ground intercept to the 

point source (Figure 1). The point source lies a distance R from the 

line-of-sight ground intercept. The vertical altitude is h, kilometers, and an 

average scattering coefficient 9g, 1/kilcmeter, assumed ccnstant, gives an 

optical depth, based on vertical elevation, of 3 h. (Optical depth (dimensionless) 

s 

based on scattering, Tg, is defined as t_ = / 6_ dz. The final program in 

s jo s 

y^pendix II includes the effect of a variable scattering coefficient, as 
a function of altitude z. ) Scattering by the atrojsphere is assumed to be 
isotropic. Except for the point source, the surface is assumed to be 
black and non-emitting. Letting I denote the intensity at the location. 


and in the direction and sense of the sensor, a solution was found in 
Reference ( 4 ) for I/I^ as a function of the dimensionless variables 
6, R/h and $gh. The solutions presented in Reference (4 ), for optical 
depths to .1 are reproduced in Figures (2) through (5). The program used 
to corpute the data plotted in these figures is included as part of the 
final progr-am, in Appendix II, in order that additional data of this type 
may be obtained if needed. Presented in dimensionless form, the results 


are valid for any values of R, h, and 6^ corresponding to the wavelengidi 
interval under consideration. 

a 





Although the small optical depths appearing in Figures (2) through (5) 
are representative of the, scattering behavior for a large range of cases 
encountered in practice , it is necessary to extend the previous results to 
larger optical depths in order to discern more general patterns of the 
effect of the point source "disturbing” emitter. Accordingly, the data of 
Figure (2) (6=0) was extended to larger optical depths and presented 
in Figure (6). Some interesting results are noted, as follows: 

(1) For small optical depths (<.l) the curves are relatively 
flat, and ,I/I q increases with increasing optical depth. 

(2) With increased optical depth to around 1, the curves show 
a quicker drop-off with R/h. 

(3) As the optical depth continues to increase beyond 1, 
however, the maximum value of I/Iq (at R/h=0) begins to 
decrease with increased optical depth. Finally, as shown 
in Figure (6), I/Iq vanishes for very large optical depths 
as it did for very small optical depths. 

(4) For larger values of R/h, the maximum value of I/Iq occurs 
with smaller values of 0gh. 

The value of Pgh for vdiich I/Iq attains its maximum value is of some 
theoreticcil interest. No atte.mpt was made in this investigation to 
precisely determine this value; however, an order-of-magnitude estimate 
indicated that for this vertical line-of-sigJ»c case the maximum wculd occur 
for Bgh<2. The cxistance of a value of 6gh for which the influence of the 
"disturbing" point is maximum is expected, and may be illustrated by the 
following example: Consider as a "sensor" a person standing on earth 

looking upward, not at the sun, cmd regard the sun as the "disturbing" 
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point source. If there were no atmosphere (i.e. , very small optical depth) 
no energy would be scattered into the "sensor” line of sight (i.e., 

I/Iq would be zero). Nex"t consider a clear atmosphere. The optical 
depth is moderate due to scattering, and 1iie "sensor" receives considerable 
visible radiation (i.e., I/I^ is greater than zero).- Finally, as the sky 
beccmes overcast, making 0h larger, a decrease in I/I^ is noted, indicating 
that "the "peak" value has been passed. 

SCATTERTO ATMOSPHERE . RESULTS FOR nNITE AREAS 

The versatili-ty of application of the results obtained for the, 
point source is evident in the extension of -the investigation to consider- 
ation of "disturbing" emitting areas of fini-te extent. Several emitting 
areas are, for conputationol purposes, replaced by point sources, each 
located within the boundary of an original area, and each point source 
assigned an "effective"' intensity erg/sec-meter^-steradian-meter. 

The effective intensity assigned to a point source is determined by the 
requirement that the point source have the same total hemispherical 
rate of energy emission in the same frequency interval as the original 
area it replaces. For a diffuse emi-tter the total hemispherical emissive potcer 
e (energy rate per area per frequency interval, erg/sec-meter^-meter) 
is related to the intensi-ty by the relation e = ifl (Reference 5). Thus 
the total rate of energy emission per frequency interval becar.es eA = ttIA, 
where A is ihs emitting area. The point source replacement cf the finite 
area is made by assigning the value obtained by normalizing the 
"effective" emitting area to unity, maintaining the saiiic hemispherical 
emissive pewer, i.e. , Igffd) = IqA, vtfiere Iq is the intensity corresponding 



to area A, This process will be exact as the size of the original emitting 
area approaches zero (or more precise.ly as A/h^ ->-0). In t}-.o canputationai 
procedure, A approaches zero' numerically by an integration tech;rLque. 

Figure (7) illustrates the method. A given area A is assumed tt.- 'on- 

unifom intensity . A four-division approximate integration of 
the contribution due to emission from that ai ea is indicated in the 
figure, in which each area segment AA^ is replaced by a point source having 
an effective intensity = Ii^^> located within the boundary of the 

corresponding area segment. Each point source is then treated by the 
methods previously developed, noting that the point sources in general 
differ in their values of R/h and (J> with reference to the sensor 
line-of- sight. The extension to * larger number of divisions for areas 
over vhich the intensity varies is clear. Numerical examples will be 
presented in a later section illustrating the methods developed above. 

Ihe above discussion eexisiders the influence of emitted energy from 
points not at the sensor line-of-siglit ground intercept. If, as is the usual 
case, the ground intercept point (i.e., the ’’target” point) is emitting . 
e;ijer^, the sensor will receive energy' directly frem that poiiit. This 
signal will be affected by the atmosphere, ancl will be reduced in intensity 
by a fraction ^ knewn as the atmospheric ’’transmittance”. In 

addition, the "target” itself may act in the role of a ’’disturbance", 
since same of its energy emitted in a direction not along the line-of-sight 
may be scattered into the sensor by the atmosphciv.. 'Ihis "disturbance’' aspect 
of the "target” may be treated as a special case of a "disturcance" for 
vAiich R=0 in Figure (1). 
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NUMEPICAt. METITODS 


A discussion of liic numerical method employed to solve the integro- 
differential equations governing the scattering problem was^ presented in 
Reference (4). Some results pertaining to the numerical method itself will 
new be presented. Evaluation of tiie multiple integrals appearing in the 
governing equations is quite tinte-consuming on the digital computer. 
Consequently, methods were sought vMcli would evaluate fixed- limit 
multiple integrals while optimizing accurary and computation time. In the 
course of the investigation a very usefxil iterative method of multiple 
integration was developed and employed to obtain some of the numerical 
<kta. This nethod is presented in ^pendix I. After development of the 
method it was determined that th extreme accuracy it provided was not 
necessary by comparison with the number of significant figures carried 
in the remainder of the conputational process, and therefore the iterative 
method does not appear in the final program in ^pendix II. 

ABSORBIM^/DglTING ATMOSPHEi^ 

Ihe investigation has thus far dealt wilh a scattering atmosphere, 
that is, one which neither absorbs nor emits radiant energy. Absorption 
and emission by the atnnosphere are of great importance in remote sensing, 
and the governing equations for a non-scattering atmosphere are 
considerably simpler than those for scattering atmospheres. In such 
an atmosphere", only those processes of emission and absorption which occur 
cn the sensor line-of-si^t need be considered, and there is no contributicn 
to ernrgy received by the sensor from any "disturbing" ematter not on 
the line-cf -sight. 
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For the non-scattoring atmosphere, Reference (&) or (6) readily 
gives the equation for the intensity received at the sensor location , 
in the direction and sense of the line-of- sight as follows: 


I = I e-^a/cos6^ _1_ 
o 4ir 



€b e 


,-(T^t)/oos8 dt 


cose 


[ 1 ] 



where Iq now represents the intensity at the ground level line-of-sight 
intercept, e^ denotes the black- body hemispherical emissive po^wr, and 

r ^ 

T«represents the optical depth based on absorption, i.e., tq= j B. dz. 

J 0 ^ 

The symbol denotes the absorption coefficient, which may in general be 
a function of altitude z. The symbols 6, z, and h have the same interpre- 
tation as in Figure (1), and t denotes a durrmy variable of integration. 

The hemispherical black-body emissive power is a function of the 
wavelength at which the intensity is conputed and the absolute temperature 
of the atmosphere. Reference (7) may be consulted for a theoretical 
discussion of this i\-nction. In order to obtain the value of the intensity I 
received by the sensor it is necessary to perform a numerical integration 
as indicated in the above equation, given 6, h, and the variation of 
the absorption coefficient and tempemture as functions of altitude. This 
numerical integration is perfonned as part of the final program in 
Appendix II. 

CONTRAST AND CONTRAST TRANSMITTANCE 

IVk) of the most importairt aspects of ;< ■ -te sensing of non-uniform 
targets are contrast and contrast transm5tt\i:vc»s. Given a "target” and an 
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adjacent' “background" j denoted by t and b, respectively, the ccn-trast at 
ground level, based on intensity at a specified waveleng-th, is defined 

to be 

C(0) = I^(0)/Ij,(0) , [2] 

idiexe the (0) denotes intensities at ground level, and the "target" may 
be assumed, without loss of generali-ty,to have the hi^er intensi'ty. Let 
I^(h) and Ij,(h) denote intensity at the sensor, at altitude h, when the 
line-of-sight passes at ground level through points t and b, respectively. 
The contrast at altitude h is given by C(h) = I^(h)/I^j(h). In general 
C(h) will be smaller than C(0) due to atmospheric participation, including 
the effects of "disturbing" emitters on the surface. The reduction in 
COTtreist is specified bj^ the con rast transmittance CT, defined as 
C(h)/C(0), and it is the contrast transmittance which measures the degree 
of pattern recognition possible. Miller, et al, (Reference 8) have 
presented some typical measured con-trast transmi'ttance values and discussed 
additional aspects of multispectral sensors for ERTS A 6 B. 

APPLICATION OF RESULTS TO REMOTE SENSING EXPERIfCNTS 

A review of ongoing and planned ERTS and SKYLAB capabilities and 
programs (References 8, 9. and others) was made, with respect to potential 
input fi.'cm the results of this investigation. A numbei' of conclusions 
and reccmrrendations result from this reviei^. Loss of contrast of a ground 
signal, may result from scattering, absorption/emission by the a'tmcspherc, 
or a combination of these effects. In some regions of the Wc'.velength 
spectrum one may be fairly certain that one effect or the other dominates 
the atmospheric attenuation (for example, attenuatior in the visible 
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range is almost entirely due to scattering). In general, ha/ever, it 
cannot be expected that the attenuation will be positively identifiable 
as to mode. In this regard it may be noted that the developient in tJiis 
investigation has assurrted that the atmosphere is either purely scattering 
or is totally non-scattering. This distinct separation was made to 
streamline the application of the results, in view of the practical 
in^ility to identify in general the mode of attenuation. Several 
ijiportant parameters can, hwever, be accurately measured or predicted. 

Ihe sophisticated 4-D atmospheric model of Chang (Reference 3) provides 
reliable tenperature and other property profiles, as well as total 
transmissivity. In addition, total transmissivity is conceptually simple 
to measure experimentally (although it is in general not possible to 
identify whether the attenuaticxi is due to absorption, scattering, or 
both). An experimental procedure for use in conjunction with a contrast 
corputational procedure might proceed as follows: Establish two 

"calibrating" emitters on the ground, emitting radiant energy of different 
intensities in the same frequency range, and positioned in the proximity 
of the areas beP^jeen which it is desired to determine contrast by remote 
sensing. These adjacent ground emiLters would be remotely sensed, and 
the transmissivity thereby determined, as the ratio of the difference 
between tiie remotely sensed intensity when. viewing, in turn, the two 
calibrating emitters to the difference between the known ground values 
of intensity’ of the emitters . This would be valid for any conbination of 
scattering/absorbing/emitting attenuation, since the atmospheric emission 
would subtract out due to the taki.ng of a difference ratio. Knowing the 
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total transmissivity, the limits on the contrast between the t\-JO emitting 
areas in question may be determined by (a) assuming the atmosphere 
to be purely scattering and tiien (b) assuming the atmosphere to be 
non-scattering. In both Cases (a) and (b) the measured or estimated 
transmissivity is used. Then the contrast and contrast transmittance 
for a partially scattering atmosphere nay be shown to lie between the two 
limits found in (a) and (b). If the predicted contrast transnittance 
range is unacceptably low, considering the kncwn or estimated intensities 
of "target’’, "background", and any "disturbing" emitters which may be 
present, the collection of extensive data could be rescheduled. 

CCMPUTER PRCX3PAM 

Hie digital coiputer program, listed in Appendix II, accepts as 
irq>ut data the gecmetry of the sensor line-of-sight with respect to one 
or more "disturbing" emitting areas whose "effective" point source 
intensities are also included as inputs. The "target" intensity is 
also read in, as is the sensor altitude, and the temperature profile 
cxarputed by the 4-D model of Reference 3. In addition, the program will 
receive as an input scattering and absorption' coefficients c>s i'unctic«is 
of altitude. The program coiputes the "band" in which the intensity at 
the sensor location may be found, by assuming the atmosphere to be, 
respectively, purely scattering and non-scattering. Hie scattering 
coefficient is not necessarily constant with altitude as was llie case with 
the preliminary results shewn in Figures (2) tlirough (5). Nor is the 
absorption coefficient necessarily constant witli altitude. 

If, as discussed in the proceeding section, scattering and absorption 
coefficients are not specified as functions of altitude, but, instead, the 
total ti\insniissivity is knavn, an avenigc value of the r>calterJ"r or 
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absorption' coefficient may be read into the conputer program, the value being 
determined from the known, transmissivity. Example calculations in the 
foliating section illustrate the coiputational procedures. 

Ihe diagram in Figure 8 summarizes the inputs and outputs of the 
ODTiputer program. 

NUMERICAL EXAMPLE 

Ihe following example illustrates the use of the results previously 
obtained, to estimate the effect of atmospheric participation upon contrast 
and contrast transmittance of upwelling radiant energy, in which the 
surface distribution of intensity is non-uniform. Figure (9) illustrates 
a typical situation of this type. It is desired to conpute the. contrast 
between emitting areaAAj^ (the "target") and areaAA2 (the "background"), in 
the presence of a participating atmosphere and a third "disturbing" emitter 
AA3. Contrast is to be carputed between t*wo adjacent points of areas 
and AA2 at the position shown as the sensor line-of-si^t ground intercept . 

For sijiplicity of calculation each emitting area is replaced by a single 
point source (rather than a number of point sources as in Figure 7) having 
an "effective" intensity based cn the intensity of the area and the size of 
the area. Three wavelengths are selected for computation, one in the visible, 
and the other two in well-knaim atmospheric "windavs". For purposes of 
ccrrparison among the ‘diree wavelengths, the respective intensities of the areas 
are taken to be the same at all wavo].engths considered. This, of course, 
would not be the case if upi-;elling radiation were due solely to the area's 
tenperaturej but, as discussed earlier, for reflecting areas the upwelling 
radiation is not necessarily an indication of the surface temperature. The 
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following table gives values used for all three wavelengths: 


1 "target” 

2 "background" 

3 "disturbance" 


2 

1 

5 


I . , erg/sec-m^ -sr-m 

lO 

. 6- X 10^** 

4 X 10^** 

- 6 X lO^** 


Average Distance 
Fran L.O.S. ground. intercept ,m. 

12 X lO^** 1 

4 X IQl** .5 

30 X lO^** 5 


(1) X = .5y 

The total atmospheric transmittance at this wavelength may be estimated 
(Reference 10) to be 0.6. Assuming this to be the only information 
available regarding atmospheric attenuation (i.e., assuming a detailed 
correlation between altitude and scattering coefficient is not available 
for input to the ccnputer program), an average scattering coefficient 
pg = .0144 is read into the progTc»n. This value of 8g is chosen to provide 
an overeill transmittance equal to 0.6, since 


Transmittance = e“'^s/cos6 _ ^-3sh/cos6 _ 


g-. 0144X25/00345° 


0.6 


It is first assumed that the atmospheric attenuation is due solely to 
scattering. When viewing the target, the magnitude of intensity at the 
sensor location is due to two contributing sources: First, the intensity 

at :'.ie sensor is equal to that at the target multiplied by the transmit Lance 
factor of 0.6, which represents the effect of energy scattered fran the 
sensor line-of-sight by the atmosphere. Second, seme energy is scattered 
into the sensor line-of-sight because of the "disturbing" emitting areas, 
represented by effective point-source emitters. There are three such 
"disturbing" emitters, counting the "target", aAj^, "background", aA 2» and 
the "disturbing" emitter AA3. Since R/h is very small for all three 


13 


i 



effective point sources, the computer program (y^pendix II) gives 
1 / 1 q = .0316 for each effective point source. (It may be noted that this 
result cannot be obtained frwi Figure 4 bec<>use the optical depth of . 36 is 
not within the range of the data presented in that figure). 

Letting denote the intensity corresponding to the sensor location 

anJ orientation when viewing the target, the above-mentioned effects are- 
combined as follows: 


I 


target 


= e-’B/=os6 , a/i^) j 

= 6 X + .0316 |^12 + 4 + 30 jlO^** 


[3] 


= 6 X lO^** X .6 + 1.4536 X 10^“ 


= 5.0536 X 10^ ^ erg/sec-m^-siMn 

When viewing the adjacent background, the energy received due to •'disturbing" 
emitters is unchanged, but the direct radiation is new received from area 
AA 2 . Il^us the intensity at the sensor, when viewing the background, is 

. Ijl«ck = ^20 [ Ileff ♦ l2eff * ^3eff J ^ 

= 4 X 10l“X .6 + 1.4536 X lO^** 

= 3.8536 X lO^** erg/sec-m^-sr-m 
The ground contrast is 

C(0) = I^(0)/Ijj(0) = Iio^Ijo = ^ 10^“)/(H X lOl**) = 1.5 

The contrast at the sensor altitude of 25 kilometers is 

C(h) = I^(h)/Ijj(h) = (5.0536 X 10^'*)/<3.8536 X lO^**) * 1.3114 
The contrast transmittance is 

CT » C(h)/C(0) a 1.3114/1.5 a 87.4% 
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If the atmosphere is nejct assumed to have a transmittance of 0.6, but 
be non-scattering, the ccmputer program evaluates intensities, contrasts, 
and contrast transmittances according to equation [1]. The results of- this 
coiputation reveal that the ccxitrast is essentially unchanged, indicating 
that the atmosphere does not emit appreciably at this wavelength. Thus 
the contrast transmittance for a non-scattering atmosphere is approximately 
100 %. 

The results are summarized as follows: at a v.»avelength of 0.5n, if 

the total atmospheric transmittance is estimated (or measured using 
calibrating emitters) to be 0.6, then the contrast transmittance for the 
above-described intensities, areas, etc. lies in the "band”, 87,4% < CT < 100%. 

X = 4y 

Chang (Reference 3) gives an atmospheric transiriittance of .86981 for 
this wavelength for Region 6 (January) for a cloudless, non-scattering 
atmosphere. For the purpose of illustrating the computation of the contrast 
transmittance "band", however, the atmospheric transmittance vail be 
assumed to be .86981 for a purely scattering atmosphere as v;ell as for 
a ncai-scattering atmosphere. If it is first assumed that the atmosphere 
is purely scattering, an average Bg = .00394 is read into the program 
to produce the desired atmospheric transmittance of .86981, as previously 
described. The program gives I/Iq = .0109 for the three "disturbing" 
effective point sources. Using equations [3] and [4], the intensities 
at the sensor when vi«-/ing, respectively, the target and the background arv 2 
confuted to be 

^1 target " ® ^ ^ [l2 + 4 + 3o] X lOl** 

= 5.7214 X 10^** erg/scc-m^-sr-m, 

I|b^ck ' ^ ^ .86981 + .0109 [l2 4 30 ] ‘X lO^** 

= 3.9814 X 10*** cr'ij/sec-m^-sr-m . 


l.'i 



The ground ccsrtrast remains at the value 1.5 for this example. The ccntrast 
at the sensor altitude of 25 kiloneters is 


C(h) = It(h)/Ijj(h) = (5.7214 X 10^**)/(3.9814 X 10^**) = 1.437 

The ccaitrast transmittance is 

CT = CCh)/C(0) = 1.437/1.5 = 95.8% 

If, next, the atmospheric transmittance of .86981 is assumed to be 

measured in a non-scattering atmosphere, the ccmputer program emplcys 

equation [1] to evaluate the intensities at the sensor location. The 

following results are obtained: 

I = 6 X lO^** X .86981 + .49604 X 10^1 

target 

^ Lack = ^ ^ .86981 + .49604 X 10^1 

As in the previous case, it is clear that the atmospheric emission (represented 

by •the second term in tfie above equations) does not greatly affect the 

signal, so that the contrast transmittance is 100% for this case. The 

ccntrast transmittance ’’band" is thus 95.8% < CT < 100%. 

(3) X= lOy 

Reference 3 gives a transmissivity of .68 for this wavelength for 
Region 5 (January) for a cloudless, non-scattering atmosphere. As in ■the 
previois example, -the transmissivi'ty of .68 will be assumed to describe 
both modes of attenuation, in order to illustrate •the coiputation of a 
contrast "transmittance "band”. As before, assuming a purely scattering 
atmosphere, C(h) = 1.347, and CT = 89.8%. For a ncn-scattering atmosphere, 
C(h) = 1.487, CT = 99.2, therefore 89.8% < CT < 99.2. 



CX)NCUISIC:IS 


The scattering, emitting, and absorbing behavior of the atmosphere is 
complex, and few gener^zations are possible within the visible or infrared 
spectral regions. Much effort has been devoted to development of atmospheric 
models to describe this bdiavior. In this connection, Reference 3, 
for example, has alreacfy been mentioned. The work of Ihrelkeld and Jordan 
(Reference 16) discusses in detail atmospheric transmissivity' on clear days. 
Threlkeld and Jordan's results indicate that water vapor and dust scattering 
are significant well into the infrared region. If detailed atmiospheric 
models are avaiilable, they may be incorporated into the program in Appendix II 
of this report. In many cases, hajever, local effects (for example, the 
discharge from a smokestack, clouds, or unusual gas concentrations), may 
require the use of calibrating emitters, as discussed in a previous section 
c this report. 

Regarding the selection of wavelength bands to employ for remote sensing, 
again no firm generali 2 ations can be made. TWo aspects should be considered, 
however. First, it is of course desirable to sense in a "winda.-?", where 
transmissivity is high. Such windows are well-loiown for standard atmospheres. 
In addition, it is necessary to ccxisider that the emissivity of the 
atmosphere varies with wavelength, and thus will produce greater interference 
with the desired signal at some wavelengths than at others. As a rough 
estimate of the wavelength variation of the atmospheric emissivd.tv, it is 
useful to c^ipute the wavelength at which the black*body emissive poi'?er 
is greatest, correspcaiding to an ayerage atmospheric temperature. It may 
be shewn (Reference 5) that XT = 5216 gives the wavelength X, microns, 
at v/hich a black body having temperature T, ®R, has the greatest emissive 


pofr/er. 


17 



ACKNOTLEDGrCNTS 


Ihe authors wish to acJqicwledge the generous support and assistance 
of NASA Contracting Officer ^presentatives, Mr. Paul Larsen and Mr. Clark 
Brown} throu^out the course of this investigation. The excellent technical 
'typing and other clerical woric of Miss Susan Shaip is also gratefully 
adoiowledged, 


18 


APPENDIX I 


EXTRAPOLATIVE METHOD OF 
NLMERICAL IffTEGRATION 



DCTRAPOLATION MCTHOD FOR MUIJTIDIMENSIOMAL QUADRATURE 


A matrix formulation of an extrapolation method for the evaluation 
of multiple integrals, applicable to any order integral, is shown to 
be concise and easily prograircned. Moreover, it is demonstrated to be 
efficient fron the standpoint of accuracy versus number of operations 
p erf ormed. The method is limited in this investigation to the evaluation 
of triple integrals with fixed limits, although it may easily be adapted 
to integrals of higher or lower order. The trapezoidal rule is enplcyed 
as the fundamental quadrature, althou^ the general multidimensional method 
can be easily developed usiivg a Simpson's or other rule. Additional 
exanples and analysis of extrapolative methods of this type may be found 
in references (11) through (15). 

Using the triple integral as the basis for the development of the 
extrapolation method, we assume that a numerical approximation is desired 
for 


I 


ba b2 bj 



fCx,y,z)dxdydz, 


*2 ®1 


where f(x,y,z) is a continuous function of the variables x, y and z in a 
volume V, and where the numbers a^^, a 2 » a^, b 2 » b^, are lower and upper 

limits, respectively, on x, y, and z in V. 

Using the trapezoidal rule, the approximation to the single integral 
is vn?itten as follows: 


I’J 



i 


I 


I,<y»z) = / f(x>y,z)dx - ^ ^ 

i=l 


y. *^) - • 


<V^i^ 2ll 

12 n2 3x2 


,a, < t < b,, [1] 


c»y»z 


wbe.re Wj^ are the trapezoidal weight ^\inctions, n is the number of intervals, 
and N = n + 1. In this equatior; the equations to follow, existence of 
derivatives to the requii^<l Cir;-.” ’In V is assumed. Formally applying this 
trapezoidal rule to a douL.v i.'iiegral, it follows that an approximation 
to the double integral my be vnc'itten as follows: 


IjCz) 5 



f(x,y,z)dxdy = > w^ / f(x^y,z)dy - 

i=l iaj. 


(b^-a^)3 p2f 

12 n2 j 3x2 


dy 

c>y»2 


[ 2 ] 


Similarly, formal application of the rule to a triple integral yields 
bg bj bj^ b3 b2 


I r 



^^2 


f (x,y,z)dxdydz = & h j^,y,z)dydz - 



i!L 

3x2 


^3 ^2 ®1 

Next, substitutions of the single integral expression are mad 2 fron 
equation (1) into equation (2), and expressicxis for the double integrals 
iq>pearing in equation (3) are replaced by using equation (2). After some 
manipulation of the expression resulting fixm this multiple substitution, 
the follov.’ing relation is obtained: 

bj bj 


dydz [3* 

C,v,z 



N N 


N 



isl j=l k=l 


w.w.Wj.f(Xi,y.,Zj^ 



?!£ 

3z2 


d>;dy 

x,y,6 




i 



7 


4 


^ 2 Lf| 

( 12)2 / 3 x 2 3 z 2 


dy - 

c»y»e 


(12)2 / 3y2 3 z2 


dx 

x,a,e 


^^3^3^^^^2~^2^^^V^1^S2 32 32 


( 12 )^ n* 


3x2 3y2 3 z2 


(bj-aj)^ 


^bi 



dxdz 


x,a,z 


( 12)2 n«. 



rr,jt 

1 3x2 Jy2 

:12 n 2 / 3 x 2 

s 

t»a,z a a 


dydzj 

t»y»z 


« . < bj 

aj < 6 < bj 


[« 


The first term on the right of equation ( 4 ) represents a trapezoidal 

rule quadrature of the triple integral, where the mrrbcr of intervals n is 

the saiTvc for the three axes. To derive an extrapolation forrxila, it 

is assuned that four separate quadratures are performed, using four 

different nunber of intervals n, viz., n^^, n2» n3» n4« The symbol 

will denote the standarxl trapezoidal approximation for the ith quadrature, 
i 

using n = n^ intervals. It is next assumed that the terns multiplying 

1 , i and ig in equation ( 4 ) do not diangc as n is varied through n2» 
n' n^ IT 

etc. With this assumption, the four different quadratures may be written 
as follavs: 
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1 




where C2> C3 denote constants. 

These equations may be written in matrix form as follot;s: 


1 




1 


1 


1 


H/i 


n. 


nf 


n. 


1 

"1 


n: 



n. 


1 


nr 


I 


^ ■ *> 
s 




C2 


In 

"3 





[ 5 ] 


[ 6 ] 


Solution of these linear equations will yield I, t^Lthin the limits of 
the validity of assuring that C^, C2 and C3 remain constant as n is 
varied. It is observed that the relevant equations for a single emd 
double-intejjral extrapolation formula may also be obtained from equation ( 6 ) . 
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For example* equation (6) produces the single-integral extrapolation 
formula if C 2 and are replaced by zeros. Thus if I is interpreted 
to be the extrapolated approx mat ion of a. single integral, its value 
is obtained from the solution of 




*4 





I 







1 

1 





1 1 




V 2J 

2 




a redundant set of equations involving I and I being disregarded. 

ng 

This is recognized as the Romberg extrapolation formula. The extension 
of the method to higher-dimensional integrals is also evident from the 
form of equation (6), in whidi the extrapolation formula for an mth 
dimensional integral is obtained by replacing C^+ 2 .’ Sn+2 • • • zeros. 

Numerical results of this method applied to a triple integral will 
now be presented. It may first be noted that, except for round“off error, 
this extrapolation method is exact for monomials of the form x^y^z*^ v;here 
a, b, and c are < 2, regardless of the choice of n^^, n 2 , n^, n^. The 
method was employed to approxirate the integral 

2 2 J2 

dxdydz 

which has the exact value. 260.002. (Thar, value and the follov5n,(» data 
are rounded off in the 6th significant figure). Hae following table 
presents the results of the computation: 





n 


I'mpezoidal Quadraturo 


Extrapolated Value (n^^, n2» n,, rij^) 


4 

277.373 

8 

264.894 

10 

263.416 

12 

262.616 

16 

261.821 

20 

261.454 

24 

261.255 


260.802 ( M , 6 , 8 , 10 ) 

260.802 ( 12 , 16 , 20 , 24 ) 


It is clear that in this case the extrapolation method is desirable. 

Hie accuracy of the standard trapezoidal quadrature degenerated due to 
round-off error, as n was increased, before the accuracy of the extrapolation 
method was approached. 



APPENDIX II 

LISTING OF COMPINTER PROGRAM 




U U U U W u u 


niMENSKirj W(35) ,C(e,?6) ,TA( ?7) ,TS(?7) ,7.<27) 

A=Cn5i(T»‘i:TA) . I- = ALTI 'UOh.KK. 

XIT = TA«>.GLT IMt-NSITY A. GROUNH, F!5r./SEC»K«*2«Sft*M 
jaPrIiACKGKCUi.P INTtUSITY AT GKC’UNCt EKG/.SEC*M«*2*SP^«.'' 


NSsTiiiMfinn OF ri STimi.i:4G fm tiers 

G=rANGIE PHI LOCATING lUL-PT' fH P I STO'O* I NG fc^.lTUR 
RP=ni STAI.CF. or IRCPTTH niSTliUBlKG r-MITTER FKOF' TARGET 
XlDlS-cFFECTlVE ItJTEAiSlTY OF. IREPT'IH DlSTURAPlG FFMTfFR«_- 


XINT=0.0 

CALL FCUROIU) 


•hRG/ScC*M**2«'SR*F‘. 


fi=26 


REAO(l» 1)A,H,XIT,X1G,*JS 

- l.f0RMAT(2F10.6;2Fl5.6,12) 

CALL OPDEPIH,fl,CA,CS,G,TA,TS,Z) 

DO . .8 . 1 = 2 , N 

rc-4=f.r.io,j j'i'r.xiM -ic'' -TAun/Ai/A 

XIf4T = XINT*FCr,«( TA( I )-TA|jn 


B. CONTINUE ... 

— 01flT = XlNT/T. 1A1S<9 

: XIA0ST=X1T/E.XP(CA/AMD1NT 

X1A0S-3=XIH/EXPICA/A.)+C1MT 

0CLTA=0.ni 


PI = 3, lAlS'i 

. RTA=S0.RT(l.-A<=O2.0l . 

_....DU 99 IRCPT=1,NS : ■ 

E=O.OOOl 

REAn(l,35)A,R0,Xin!S 
35 FnnMAT(2»U0.6,El5.6) 

0=R0/H - ...... 

xxx=o.o :.... 

DPH1=2.«3.1A159/X.\ 

DXMll=( l.-UELTA)/XN 
no TOO It=2,N,2 


..XIMT^O.O • 

COMl = E/ A«l‘2.0*0<=*2.0/E-2.0«04'RTA*C0S(B)/A 

CCN2 = 2 -O^F..-.. 

r.0N3 = ?.0*E^<RTA 
CTE=T$(JE) 

_.XKU=DELTA ' 

6..PHl = CUO 

RTXF*U=SCRT(l.-XKU<'#2.r»)__ 

VAH3*F/l XHU**2.0) 

VArXsP.OoOORT XMU 

— VAR7»CON3«RTXNU .... . 


JJ=l . 

A VARl-VA.RbfCCSIPHn/XyiU 

..VARA=VAR70CCS(P-PHn/lA#XMUl ... 

YAR'>=a)N?-VARA 
1 = 2 


2 J=l-l 

VAR2»l.-7( J) /(HtE) 

IF(CrC>TSUn50,50,52 

_ . _50..V'Rl--VAR6^CUStPHl)/XMU ... . — 

VA!<5=C.IN?+VARA 
02 CrjU lNIIr 

TLRM2»t*(CnNl-VAR2<-( VARl“VAR3*VAR2+VAR5) ) 
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R = TS(J)'!‘H«S0'^T(rrRf'.2>/Z(J) 

xAr>=M<S(Cni-is(j )) 

TLST^R4XAR/X.VU 

ir ( TrsT. .CT. rj.fuGO tc 75 

. . FCN=l.O/(£XFMtUXAli/XrU)*XI>lU) . . .... . 

IFU.EQ.K)M = N ' 

X1.NT = XI,N1 +t)X«U«OPHI.’fL( T5(.M)-T.S.CJ1 ).*FXN 

75 COMTINl'F 
I = I+? 

1F(I-N)2,2,3 . 

PHI = PHI+DPHI 1 

JF.( JJ-'J)A,A,3 

5 K=K+1 ' 

XyU=XMU+CXVU 

IF(K-.\)5,6,7 . 1 

7. .CONT.INUE 

XIEX=XINT ' 

T£RMl = (.(RTA*CUS<e)/A-D/E).«lRTA*CQS(B)./ArD/.e) + lRtA*SI.NlB)yA)**2.tl 

1.0) OX'0 . 5 

G=1.0/(EXP(CTE«TcR.Ml ) )+X I CX/ .0*P I) 

. ARG=G/lAtEXPnCS-CTE)/A)*A.*Pn : 

IF( lE.EC.iM) I r=lE :. .'. - ■ 

XXX = XXXt ARG^'(TSnT).-r${J.E).) ^ ^ 

IF( JE-l ) 55,55t-56 
55 F=E+,0H12 H3 H 

.... .56...E = E + .0b3 3333i3 ...! 

._...57 .CCNTIN'Ut 

JPO__CnNTINUii ' 

SCAT = SCAT + XXX<=XiniS 
<jq cofni JiJE 

XISCAT = XIT/CXP(CS/A)*SCAT . : 

XISCAO=XIB/tXP(CS/A J +SCAT 

....GCO.M = X1T/XIO .’. 

HCO'JS=X) SCAT/XISCA9 1 

HCONA=XI A9ST/XIA9SB ’ ‘ 

CTRANG = HCO'iS/GCON - 

CTRANA=HCUNA/GCGU 

401 FllRMAT(60X,' INTENSITY*, //) 

402 FORMAT! OX, 'NON SCATTERING ATMOSPHERE T85 SC ATTER INC ATMCSPHERE*, 

"aFi FORfUT ( T uF tFrF^* , T43, • PACK^nUNO ' , T35 ,* TARGET •, TGA. *BACKGROUNO* 
It//) 

404 FCRMATIGIX, ’CCNTRAST*,//) ...... 

.405 . FORMAT I 54Xi ’CONTRAST TRANSM I TTANCE ' t // ) 

WRlTE{3|40n ; 

-KR I T E ( 3 , 4 02 ) - - 

WRITE! .1,403) 

WRITE! 3,4 0ft)Xl A0ST,XIAPSI3,XISCAT,XISCAR 

WRITE! 3,404 ) 

WRITE! 3, 402) ' 

WRITE! 1,407)HCONA,HCONS 

WRITE! 3,405) . 

WRITE! 3,402) 

WRITE! 3,407 )CTRANA,CTKANS 

405 r0RMAr!T4,El5.6,T4l,C-l5.6.T8l,Fl5.0,T97,El5.6,//) - 

407 rURMAII T25,E15.6,T89,E15.6//) 

CALL EXIT 
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SUR’tnuriNE OPOrP(H,.N,CA,CStO,T/V,TS,7 ) 
CIMPMSICN 0(6t26) f TA(?7) ,TS(27) ,Z(27) 
07 = H/(\'-l ) 

Z(1)=0.0 

TAC I ) = o.n 

Tsn )=o.o 

CO ‘)9 1=2, N . .. 

..z( n=z( j )*f)z , .. 

_ .. FA( I J=TA( J)*'3£TA/\( J, C)*CZ . 

99 TS( n = TS( J)*RHrAS(J,C)«nZ 

Z(1)=0. 00001 

TS( I ) = 0 .O-JOO 1 . 

CA=TA(26) 

,.CS=TS(26-) 

..i<£.TUnrJ 

tNf) 


FUN'CTION fcR(U,l) 
niMENSION 0(t*,26l 
XL=WAVf:LeNGTH, PETERS 

_IK-=D(.T, I J 

XL=A.C E-6 
1=6.62'. E-27 

— C=2.99« E 8 - . 

R=l.38 fc-l6 

Efi=(2.0«3.lA159’SH«C*«2.0)/{XL««5.C*(£XP( (inC ) / ( XL’S'B’t'TK 



RET UR*'.' 

EK'D 


FUNCTION RETAA( 1 ,0) 

CIP.tNSIUN r.M6,26I 

PETAA = AP.S0RP1 HJN COEFFICIENT, l/KP 
UETAA=0.C03?A 

RETURN 
ENO 


FUNCTION RETAS(I,D» 

DIMENSION 0(8,26) 

eFTAS=SCATTEP.i:iG COEFFICIENT, l/KP 

__I5ETA3=0.C06 ; 

RETURN 

END 
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SUPP.OUTIIMC- FOU-^DCO) 

C cor.PUTE-cu'^vES OF ‘atmosphere properties 

DIMENSION CRT (7) ,01 8,?6),YY( 11), ZZ (60) 

GROA.2 . 

Rl=^A6.3e 

204 i\C = ‘) 

- 2C0 FORMAT (12) . . 

-Nl=l 

....\2*NQ 

,C READ IM ZZ ROW ■ ' • ! 

205 DO 201 J=l,8 
Nl=5*J-4 

- . - _ N4=N3 + 4 . ■ . . . : 

201. READ (1,600, ERR=9<59UZZ(U ,L=M3,N4.) 

6C0. FORMAT (6X,5L14.7). ' 

-601 FORMAT (2 12, 2X, 2El4 . 7, 3E 1 1._4 , 1 6, IX, 31 2J 

N3=Ni+5 

N4=N3*1 

READ! 1,6G1,EKR=999) IREG , MONTH , I ZZ ( L) , L=K3 , N4 ) , Y7, CRT (6 ) , 

1 G R T ( 7 ) , M P M , M G R , N L , N U 

— CRT ( 5 ) =NGR 

GRT ( 5 ) = GRT 1 5 ) V ICC0,.C - 

PiV.M=Nn,v, 

Pf*M=ALUG(PMM/ 100.0) 

no 540 J =2, 7 ^ 

- IF(J>7) 1113,1111,1111 : 

... nil YY(1)=Y7 

CO .TO . U12 

1113 DO 202 K=i,NQ 
N3=NC*( J-2) 

2C2 YY(K) = ZZ(N3 + U) . 

.C.COSPUTES 26. PTS IN .J-TH. CCLUMM WITH ..JH( J) .CCEF.S YY.(L ) FCR. GRT (J ) , D (J, K 

-...11 12. .CO 300 Kri,26 : 

HK=K—1 _ _ 

I F ( ( ( .1- 7 )~* ( J-67 ) . HE ,~0 ) GO ft)’ 1 7"o ! 

IF(K-7) 170,020,521 

-.-.920 ir(J.E0.7) LM7,6) = SCRT(0(7,6)) . -. 

921 D( J,K)=0( J,6)/EXP(GP.T( J)«(XK-5.0) ) 

- _...IF(J-6) 171,300,171 :. 

17-0_- 1 r ( X K . t C . 0 , 0 ) ..GC TO-. 1 1 1 

C( J,K)=0.0 
NX=NQ 

1FU.E0.6) UX=6. 

IKJ.E0.7) NX=1 

DO... eo L = 1,NX 

R0_C( J,K) = f)( J,K )+YY(L )«XK«*(L-: 1 

CO TO 652 
111 D(J,l)=YY(l) 

.--.652 IF (J.EC.3) GO TO 653 

IF1J.LT.5) GC TO 654 

D( J,K) = D( J,K)/EXP(GRT( J)*XK) 

654. IF. (J.EQ.6) 00-TC..3CC : 

171 D(J,K)*nU,K )**2 

GO TO 100 . 

653 C(J,K)»1.0/D( J,K) 

300 CONTI.NUE 

IFIJ.NE.3) GO TO 540 

-699. CC 655 1 = 1,26 

D(l,I)sO.O 

655 C(n,I)»0.0 

DO 656 K«l,26 

XKaK-l. . 
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1 

i ' . 

if , 

1 '' 

658 CO 657 L = l, NX Z Z1 .,ZI I ’ 


657 f)n»K)-()(l,KH-YY(L)<>' XK ♦« L/XH 
Cn ,K)=fX-’(PNK-Gl?’^P( 1,K) ) 

656 CCn.KjsOn.KJfRI/DO.K) .. 

.. . 560 CCMINUE 

6 22 C(7,6) = D.(7,6)«*2 

C Ef.0..C0^*5UTL. . . 

WRITE (3, 9606 )MnNTH,0( 1, 1 ) ,D( 1 , 1 ) , Il^Cr, 

9696 FORMAT! 1HI,6U, 36HEOUR-CI..VcNSIONU A TNCSPHcR I C MODEL / / 1 7X , 7HM0NTH 

- - - — l=»13,16X,10HLATITUr;E = , F 3.2 ♦ 3X , 1 1 HONG I TUCE =»F3.2,2H C, 12 X, 6 HREG 1 

2CN, I3//1 ?X,5HL£VEL,3X,8HPRESSURE,2X,>21)VARlANCE,2X,aHTEMP (K) ,2X,8H 

3VARIA.NCE,1X, lOH' WATE.R , lOH. VAR lANCE lOH ... DENSI TY ,9H VARIANCE/ 

t:0 766 K = l,26 ' 

N=K-1 . 

.766 WRITCIJ, 9697) N » ( 0 ( L , K ) , L = 1 , 6 ) » ( D ( M , K ) , M=6, 8 ) ♦ 0 ( 5, K ) . 

9697 FMRMAT(16X,I2»3H.KM,2X,8F.1U.2) 

999 RETURN 

E-ND 
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FIGURE 7 



Establishment of Emission Point Sources in an Area 
o£ Non-Uniform Intensity 


(a) Emitting surface area>A, having non-uniform intensity. A ® ^^i 

(b) Replacement of emitting surface area by point sources, “ ^i^^i 
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INPirr DATA 


Sensor Altitude, h, (Figu^ 1). V;aveleng:th, X. 

Sensor linc-of*sir,ht inclination v;itli vertical, 6. 

4-D Atmospheric mcxiel, based cn sensor location and time of year' 

3g,6^, scattering and absorption coefficients, respectively, as functiais 
of altitude z 

Number of "disturbing” emitting sources, NS 

distance of ith disturbing source from line of si^t ground intercept 
rotational angle locating ith disturt)ing source 
Effective point source intensity of ith "disturbing" emitting area 
"Target" and Backgixxind" surface intensities, at line of sight ground 
intercept 


oiTTPirr 


Intensity at sensor location, direction, and sense, assuming a non -scattering 
atmosphere • 

Intensity at sensor location, direction, and sense, assuming a purely scattering 
atmosphere - 

Cmtrast and Conti\ist Transmittance between "Target" and "Background", 
assuming attenuation mode to be (i) purely scattering ai^d (ii) nen- 
scattering. 


HGURE 8 

INPLTT />ND OUTPUT DATA DOR COMPUTER PROGRAM 
LISTED IN APPENDIX II 
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GEOKin’RY or DWIPI^: IlJUSmTING 
COMPUT/VnON OF C0MJU*VST llWiSMIOTANa: 
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